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The  Zenith-angle  distribution  of  stopping  muons  at  sea  level,  and 
the  response  of  a stopping  muon  cosmic-ray  detector  were  investigated  in 
this  research  re nor t. 


In  order  to  determine  the  angular  response  of  a wide-angle  stopping 
muon  telescope  ■ smaller  detectors  were  placed  in  coincidence  with  it  ^ 'pro'-'W 
vidiA^t*  correlation  of  count  rate  with  zenith-angle  of  particle  arrival. 

Usinp  numerical  integration  techniques,  computer  solutions  were  employed  to 
predict  the  zenith-angle-selactive  count  rates  for  various  descriptions  of 
the  incident  flux.  - A least-squares  analysis  yielded  a cosine  powers#  . . 

n * !»,7  t 0.‘j  in  the  flux  expression  1(0,  * I()  cosn  0. 

* -•  Employing  this  flux  expression  in  predictive  computer  programs 
provided  an  accurate  model  of  the  wide  angle  telescope  response  from  which 

the  telescope  half-angle  was  determined.  A—  Considerations  of  overburden 
effects  and  data  corrections  are  discussed,. 
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ABSTRACT 


In  order  to  determine  the  angular  response  <<t  <■  wnie-unqh 
stoppinq-nuon  telescope  , snallet  deteetnrs  were  plaod  r 
with  it,  providing  a correlation  of  cxxint  rate  with  zenith  angle  >f 
uirticle  arrival,  'ising  nuneriral  integration  technkjues,  or<g*itr-! 
solutions  were-  (ployed  to  predict  the  zemth-angle-'selective  orxmt 
rates  for  various  descriptions  ot  Uk^  incident  flux.  A least -squares 
analysis  yielded  a cessine  f*iwer  of  n •-  4.7  4 0.5  in  the  flux  expression 
I > I cosr  ••.  {ploying  this  flux  expression  in  pn<lictive 
('otnfjuter  oroqrams  provided  an  accurate  model  of  the  wide  anqle  telosoop' 
response  tm  which  the  telescope  half-angle  was  determined.  Considera- 
tions of  vertxirden  effects  and  data  corrections  are  discussed. 
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A.  Nature  of  the  Cosmic  Radiation  Beir.  J Investigated 

Ever  since  the  discovery  of  cosnic  radiation  at  the  beginning  ot 
this  century,  men  have  attainted  to  answer  several  questions  concerning 
it.  The*  significant  problems  in  cosnuc-ray  astronomy  require  the*  dis- 
ciplines of  astronomy,  astrophysics,  particle  physics,  and  even  meteor- 
ology, in  order  to  answer  than.  Such  questions  as  those*  are  asked: 

Where  do  cosmic  rays  originate? 

What  mechanisms  produce  and  accelerate  them? 

What  happens  to  cosmic  rays  in  their  transit  of  deep  space? 

How  do  they  interact  with  the  earth's  atmosphere? 

Of  what  particles  are  they  oerposed? 

How  may  cosmic  rays  be  detected  and  counted  accurately? 

If  these  questions  are  examined  in  reverse  order,  it  becomes 
apparent  that  a good  deal  of  extrapolation  is  required  to  utilize  data 
collected  at  the  earth's  surface  if  events  in  deep  space  are  to  be 
analyzed.  This  extrapolation  is  oonpilicated  by  the  interactions  of 
cosmic  rays  with  the  earth's  atmosphere,  and  with  the  geomagnetic  and 
heliocnagnetic  fields. 

A diagrammatic  description  of  cosmic  radiation  is  shewn  in  Figure  1. 
Measurements  taken  directly  in  the  primary  region  indicate  that  the 
primaries  are  basically  isotropic  above  the  horizon  O,  and  are  com- 
posed mainly  of  highly  energetic  nucle- , especially  protons.  A rtore 
conplete  description  of  the  primary  atomic  abundances  is  contained  in 
Table  I^?).  In  addition  to  these  fast  nuclei,  the  primary  radiation 
contains  sane  neutrinos  and  energetic  gairma  rays. 

r 
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Atmospheric  Nucleus 


Figure  1.  Examples  of  cosmic  ray  interactions  with  the 

atmosphere. 


Of  the  secondary  commie  rays  surviving  at  the  earth's  surface, 
by  far  the  most  plentiful  constituents  are  photons,  electrons,  and 
auons.  the  photons  and  electrons  emprise  the  "soft"  cmtrment, 
which  is  separable  f ran  the  "hard''  component  by  shielding  with  approx i- 
irately  six  inches  of  lead.  The  hard  ocaponent  is  almost  exclusively 
muons,  with  protons  and  pions  comprising  only  one- ha  if  percent  of  the 
hard  radiation (''“K 

The  rtf  and  properties  of  the  muon  are  rather  well-knrwn.  It 
is  an  unstable  particle  with  a half-life  of  about  1.6  microseconds , and 
has  a mass  of  about  106.6  MeV.  It  is  bora  in  the  decay  of  two  unstable 
fra<jnents  of  a nuclear  explosion  in  the  upper  atmosphere,  caused  by  the 
impact  of  a primary  particle  upon  a nucleus  of  the  atmosphere.  These 
two  fra^ients,  the  charged  pion  and  the  kaon,  decay  as  follows ( 1 : 

it*  -*  yf  + v 

with  a half-life  of  26  nanoseconds,  and 

K-  •*  y+  + v 

with  a half-life  of  about  12  nanoseconds.  K-  and  itt  represent  charged 
kaons  and  pions,  respectively,  the  parent  particles  of  virtually  all 
ivuons  in  cosmic  radiation^,7'P’^^.  The  pion  decay  schane  is 
virtually  one  hundred  percent  probable,  and  the  kaon  scheme,  about 
sixty  percent  * 5 ) . 

The  majority  of  muon  production  takes  place  in  the  region  of  the 
atmosphere  ten  to  fifteen  kilometers  above  the  earth's  surface *P* * 76) . 
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The  ext  1 cmeiy  short-lived  pions  and  kaons  decay  in  flight,  but 
miens  have  sane  probability  of  reaching  the-  earth’s  surface,  aided 
greatly  by  the  relativistic  time  dilation  for  fast  particles.  The 
survival  probability  of  the  muon  is  defjondent  not.  only  on  this  deca> 
probability,  but  also  on  the  amount  of  energy  it  has,  which  is 
decreased  try  10m 7at ion  loss  m its  transit  through  the*  atmosphere. 

Six  hi  Id  tin*  muon  be  brought  to  rest  before  it  decays,  it  will 
remain  virtually  noninteracting  until  its  decay.  Whether  this  event 
occurs  in  air,  rock,  or  inside  a detector,  the  decay  follows  the 
reaction 


* * 
t;  ♦ e“  + v 

with  a half-life  of  about  1.6  microseconds.  As  shall  be  seen  later, 
this  decay  can  be  useful  in  identifying  a particle  as  a muon.  Figure 
2 diagrams  the  flight  of  two  muons,  one  of  which  has  so  much  energy 
that  it  penetrates  deep  into  the  earth’s  crust,  and  the  other  having 
just  enough  energy  to  penetrate  and  stop  within  a "detector"  at  the 
earth's  surface. 

The  stopfling  muon  is  the  particle  with  which  this  experiment  is 
uniquely  concerned. 

B.  Detection  of  Cosmic  Rays  in  General,  and  of  Stopping  Muons 

The  first  detections  of  cosmic  rays  were  made  with  electroscopes, 
which  were  found  to  lose  charge  anomalously  to  the  atmosphere ' 6 . 
The  identification  of  the  cause  of  charge  loss  as  ionizing  radiation 
was  first  accomplished  with  cloud  chambers;  later,  photographic 
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emulsions,  spark  charters , and  bubble  chanters  were  used  to  nake 
visible  then  fl lqhts  and  their  interactions  with  matter.  Once  they 
were  positively  identified  as  charged  particles  of  various  kinds,  the 
detection  problem  became  on  of  detecting  them  efficiently,  cheaply, 
and  in  large  nutters  without  excessive  random  error.  'The  first 
detector  used  for  this  purpose  was  the  well-known  Geiger  counter. 

Another  detection  method,  the  only  one  enployed  in  this  experiment, 
employs  scinti1 lating  materials.  The  same  ionization  which  triggers  a 
Geiger  tube  causes  a scintillator  to  emit  a snail  flash  of  liqht,  the 
anplitude  of  the  light  pulse  beinq  proportional  to  the  amount,  of  energy 
lost  by  the  particle  through  ionization.  This  light  pulse  is  detected 
by  a photomultiplier  tube,  whose  output  is  an  electrical  pulse  caused 
by  the  photoelectric  effect  with  subsequent  anplif ication.  Scintil- 
lating materials  are  of  many  types,  but  the  scint ilia tor  used 
exclusively  in  this  experiment  is  a transparent  plast lc,  much  1 ike 
plexiglas,  which  has  been  doped  with  a chemical  scintillator. 

Electronic  coincidence  circuitry  provides  the  means  whereby' 
different  combinations  of  detectors  may  be  used  in  cosmic  rav  detection 
ami  selectivity.  The  principle  of  coincidence  logic  is  illustrated  in 
Figure  3,  a diagram  of  a hypothetical  cosnic-ray  telescope.  Pj,  P^, 

P3  and  P4  could  be  trays  of  Geiger  tubes,  or  they  could  be  sheets  of 
scintillating  plastic  with  attached  photomultiplier  tubes.  A coin- 
cidence between  any  two  of  them  could  indicate  that  a cosmic  ray  liad 
transited  them,  or  it  could  have  been  a randomly-generated  coincidence, 
called  an  accidental.  T V<d  separate  cosnic  rays  arnvinu  simultaneously, 
or  noise  pulses  appearing  as  cosmic  rays,  can  cause  accidentals;  the 
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note  coincidences  that  arc  required  for  the  event  signature,  how- 
ever, the  loss  likely  such  accidentals  should  be,  since  the  upper 

> 

bound  on  accidental  counts  is  given 

W = Trn-lRlR2  ...  Rh 

whore  Tr  is  the  resolving  time  of  the  coincidence  circuitry  (typically 
a fraction  of  a microsecond) , and  the  Rj's  are  the  count  rates  lor 
each  detector  by  itself.  The  event  P2  P3  P4  has  very  little 
probability  of  being  an  accidental,  and  would  be  interpreted  as  a 
through  cosmic  ray. 

Anticoincidence  measurements  can  be  equally  valuable.  The  event 
P|  P2  P3  P4,  where  the  bar  indicates  the  absence  of  a coinciding  pulse, 
signifies  a cosmic  rav  which  has  transited  the  top  three  plates,  but 
which  stopped  before  reaching  the  lowest  plate.  S could  be  a spark 
chamber,  a cloud  chamber,  a mass  of  scintillating  plastic  tor  analyzing 
muon  decay,  or  simply  a mass  of  lead.  L,  a slab  of  lead,  is  often 
used  to  eliminate  the  soft  component  of  radiation  from  a telescope  of 
this  sort. 

TVo  properties  of  muons  are  useful  in  detecting  them  uniquely  in 
the  cosmic  ray  flux.  First,  all  muons  interact  only  weakly  with  nuclei, 
and  thus  can  transmit  large  thicknesses  of  shielding  without  being 
stopped.  As  a natter  of  fact,  much  ooanic-ray  research  is  carried  out 
deep  underground,  in  order  to  insure  the  purity  of  the  muon  component . 
Also,  the  muon  decay  into  an  electron  nay  be  employed  to  give  the  muon 
a unique  signature.  T'io  latter  method  is  of  importance  in  this 
experiment..  When  a muon  is  stopped  by  ionization  loss  within  a nass 
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of  scintillator,  two  pulses  will  be  observed  by  the  viewing  photo- 
tubes: the  first,  pulse  from  the  muon  itself,  the  second  pulse  from 

its  decay  electron.  As  an  added  benefit,  the  time  differences  between 
a large  number  of  these  pulse  pairs  should  correspond  to  the  muon 
decay  curve,  i.e.,  should  represent  a radio-active  decay  with  a half- 
life  of  1.6  microseconds . 

Thus  one  signature  of  a stopping  muon  could  be  the  observation 
of  a stopping  particle  via  coincidence  and  anticoincidence  techniques, 
followed  by  a delayed  pulse  seen  only  in  the  stopping  region. 

C.  Nature  and  Purpose  of  the  Experiment 

A wide-angle  stopping-muon  telescope  has  been  maintained  under 
the  auspices  of  the  Physics  Department  of  the  United  States  Naval 
Academy*  and  as  of  this  writing  has  been  collecting  data  for  approxi- 
mately two  and  one-half  years.  Figure  4 indicates  the  major  components 
of  the  teiesoope.  It  consists  of  a coincidence  scintillator  plate,  a 
scintillating  stop  tank,  and  an  anticoincidence  plate  to  fulfill  the 
stopping  signature.  The  purposes  of  this  experiment  were  to  determine 
the  angular  response  of  this  detector,  and  to  determine  an  expression 
for  the  stopping-muon  flux  within  and  outside  the  surrounding  building. 

Knowing  the  detector  response  as  a function  of  zenith  angle  would 
provide  knowledge  of  the  number  of  particles  arriving  from  the  north 
and  south  of  zenith,  a measure  of  the  sentient  of  the  celestial  sphere 
actually  being  observed.  It  would,  in  addition,  be  an  exact  measure 
of  the  expected  time-width  of  a narrow  unfocused  source  on  the 
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Figure  h.  Major  r\ 'ntillator  cor.ponor.tr  of  stoop » ng-muon 

telescope. 


celestial  sphere,  since  t!ie  telesocxie  is  rotated  by  the  earth  in 
a west-oast  direction.  Knowing  the  actual  distribution  of  stowjinu 
muons  at  sea  level  wrxild  be  an  excellent  test  of  the  weaknesses  <*• 
various  muon-production  models  that  have  been  develofied  elsewheie. 

The  measured  response  of  the  detector  at  various  zenith  anti 
azimuth  angles,  obtained  by  using  narrow-anqlo  detectors  as  arimt  ional 
coincidence  units  to  define  the  arrival  directions  of  *he  muons,  was 
to  be  compared  to  ocr^juter  programs  which  predicted  this  response. 

Sine"  different  predictions  of  the  measured  response  would  be  generated 
for  different  assumed  fluxes,  curve-fitting  techniques  were  employed 
to  determine  the  best  descriptions,  both  of  the  flux  and  of  the 
detector  response. 

Secondarily,  the  stoppinq-muon  distribution  obtained  was  to  be 
compared  with  differential  fluxes  measured  by  other  autliors.  This  was 
to  give  an  indication  of  detector  efficiency,  of  tiie  correctness  of  t in- 
flux egression,  ami  of  the  effects  of  the  building  upon  the  rruvasur*- 


ments. 
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II.  APPARATUS 

A.  Wide-Angle  Stopping-Muon  Telescope 

Figure  5 shows  a more  detailed  schanatic  via*  of  the  wide- 
angle  stopping-muon  telescope,  the  angular  response  of  which  was  to 
be  determined.  Each  of  the  48"  x 48"  scintillator  plates  was  viewed 
by  four  RCA  56  AVP  2"  photomultiplier  tubes,  operated  independently 
of  one  another.  The  light-tightness  of  these  plates  was  maintained  by 
aluminum  sheets  on  their  top  and  bottom  faces,  and  black  polyethylene 
around  the  edges.  The  nuan  stop  tank  consisted  of  fifteen  1"  thick 
scintillator  plates,  each  16  3/4"  x 16  3/4",  stacked  sideways,  with 
one  3/4"  thick  light  pipe  on  each  end.  The  light-tight  box  surrounding 
the  scintillator  was  constructed  of  black-painted  3/4"  plywood.  The 
stop  tank  was  viewed  by  two  Amperex  XP  1040  5"  photomultiplier  tubes 
operated  in  coincidence. 

The  signature  required  to  classify  an  event  as  a stopping  muon 
was  a TSB  coincidence,  followed  by  an  STB  coincidence  between  0.3  usee 
and  5.0  nsec  later.  The  TSB  pulse  was  also  used  as  the  start  pulse 
for  an  LRS  time-to-pulse  amplitude  converter,  and  the  SIB  as  the  stop 
pulse.  The  time-to-amplitude  converter  provided  time-scaled  output 
pulses  to  a Kicksort  multichannel  analyzer.  The  distribution  of  pulse 
pairs  with  respect  to  time  interval  was  expected  to  correspond  to  the 
muon  decay  distribution.  The  multichannel  analyzer  accumulation  was 
also  used  as  a check  on  the  ability  of  the  system  to  separate  true 
muon  counts  from  noise-generated  counts. 

The  2”  photomultiplier  tubes  were  operated  at  2200  V.  The  5" 
tubes  were  operated  at  2300  V for  part  of  the  experiment,  and  at  2200  V 
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thereafter  (see  Section  TV  E.)  The  top  and  bottom  plates'  tubes 
were  supplied  from  a high-voltage  supply  separate  from  that  for 
the  stop-tank  tubes,  effecting  electrical  isolation.  Appendix  A 
contains  the  complete  block-diagram  circuitry  for  the  stopping-muon 
counting  system. 

B.  Zenith-Angle  Detectors  and  Mount 

Figure  6 shows  a disassembled  view  of  one  of  six  zenith-angle 
detectors,  three  of  khich  were  used  in  the  measurement  of  the  response 
of  the  wide-angle  detector.  The  two  scintillator  discs  were  separated 
from  the  photomultiplier  tube  by  a non-scintillating  plastic  disc  in 
order  to  reduce  ringing  in  the  tube  from  large-anplitude  light  flashes. 

The  scintillates:  was  wrapped  in  aluninum  foil  to  enhance  light  collection 
efficiency.  The  single  Amperex  XP  1000  2"  photomultiplier  was  operated 
at  1600  V,  the  power  supply  being  separate  fran  the  stop  tank  supply 
and  the  top  and  bottom  plate  supply. 

The  basic  probes  themselves  had  been  constructed  previously,  for 
use  in  another  experiment.  When  this  experiment  was  ccrmenced,  each 
phototube  was  optically  recoupled  to  the  light  pipe  using  Dow-Oorning 
20-057  Optical  Coupling  Ocmpound.  The  masonite-aluminum  joints  were 
sealed  with  electrical  tape,  and  the  probes  were  tested  for  light- 
tiqhtness  and  correct  amplitude-’-esponse  curve.  A wooden  "key", 
measuring  4"  x 8"  x 3/4",  was  added  to  the  scintillator  end  of  each 
probe,  for  mating  with  the  radial  zenith-angle  mount. 

The  zenith-angle  mount  was  designed  and  built  subject  to  the 
following  considerations:  that  the  detectors  it  supported  should  be 
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equidistant  from  the  stop  tank  center;  that  they  should  be  able  to 
be  positioned  anywhere  up  to  the  limits  of  the  muon  telescope's  viewing 
angle;  and  that  measurements  should  be  available  in  the  north-south 
plane,  the  east-west  plane,  and  selected  diagonal  planes.  The  main 
structural  elements  of  the  mount  are  four  plywood  arcs,  of.  radius 
6Q.0",  which  are  fastened  together  in  pairs  at  maple  slot  boards  every 
ten  degrees  of  arc.  These  maple  slot  boards  also  support  the  zenith- 
angle  detectors.  The  two  arc  pairs  are  slotted  in  order  to  mate  into 
four  perpendicular  arms.  A fifth  arm,  also  of  radius  60"  but  of  only 
half  the  length,  was  designed  for  positioning  in  planes  along  a diagonal 
azimuth.  Figure  7 shows  an  orthogonal  vieu  of  one  mount  arc,  and 
Figure  8 is  a perspective  view  of  the  assembled  mount  with  diagonal  arm. 

The  original  design  of  the  mount  called  for  resting  it  directly 
on  the  top  coincidence  plate  of  the  muon  telescope.  However,  in  order 
to  clear  the  phototubes  on  the  edges  of  the  top  plate,  part  of  the 
bottom  of  each  arm  was  removed,  causing  the  depth  of  recess  to  be  only 
6",  and  the  assembled  mount  was  rested  on  four  canent  blocks  turned 
edgewise,  placing  the  recess  approximately  7.5"  above  the  top  plate. 

Figure  9 is  a schematic  representation  of  the  entire  detection 
syston.  The  automatic  data  readout  system  printed  out  five  separate 
counts  every  thirty  minutes,  namely; 

- TB,  the  through  coanic  ray  rate; 

- TSfi  * STB,  the  stopping  muon  rate; 

- Z4TSB  * SfS,  Z9TSB  * STB,  and  Z15TSB  * STB,  the 

three  zenith-angle- selective  stopping  muon  rates.  The  asterisk  repre- 
sents the  delayed  coincidence  requirement.  The  subscripts  4,  9,  and  15 


Limits  of 
top  plate 

Cement  block 
for  support 


Clot  for  m-itlnr:  v.ith 
perpendicular  irm 
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are  merely  labels  which  identify  particular  zenith -angle  detectors 
and  their  associated  electronics. 

The  dimensions  of  the  detector  system  are  listed  in  Table  2. 

The  predictive  ccnvuter  programs  explained  in  Sections  III  A and 
III  B are  based  on  these  dimensions.  One  of  the  first  discoveries 
made  during  this  experiment  was  the  offset  of  the  stop  tank  of  0.'" 
from  the  center  of  the  detection  systan  defined  by  the  top  and  tut  ten 
plates.  'This  turns  out  to  have  a significant  effect  on  the  count  rate 
of  a zenit h-angle  detector  at  the  40*  N and  40°  S positions. 

C.  Detectors  used  in  Efficiency  Measurements 

TVo  small  detectors  were  erployod  to  measure  the  efficiency  ot 
various  system  oon*»nents.  One  was  a 4"  x 4"  x 3/4"  single  scinti- 
llator, v«%kx1  on  its  face  by  an  XP  1000  phototube.  The  otlier  con- 
sisted of  two  stacked  6"  x 6"  x 3/4"  scintillator  plates,  viewed  on 
edge  by  one  XP  1000  phototube.  Each  tube  was  separated  from  the 
scintillator  by  a 1/4"  thick  non-scintillating  Lucite  liqht  pit**. 
Light-tightness  was  provided  by  electrical  tape  and  black  polyethylene 
sheet.  Springs  were  used  on  the  6"  detector  to  rraintain  pressure  on 
the  optical  }oint,  aid  its  scintillating  plates  were  wrapped  in 
aluminum  foil  to  enhance  light  collection  efficiency. 

• 


III. 


PREDICTION  OF  STTWIMV-MUON  TCLESCnpE  RESPONSE 


A.  The  Prediction  Preplan 

Orve  may  mkc  an  analog  between  the  problan  of  predict  ing  a 
detector's  resuonse  to  an  unknown  flux,  and  a boy  shoot inq  marbles 
into  a rink.  The  unknown  quantities  are  the  number  and  rate  of  marbLes 
which  the  boy  shoots  in  each  direction,  .rnd  the  nunber  that  arrive  m 
the  rink.  One  oould  place  cups  within  the  rink  for  a uniform  length 
of  time,  retrieve  the*".,  and  count,  the  marbles  within  than,  obtaining  a 
rate  distribution  over  the  area  of  the  n.x.  This  data  oould  be  used 
in  summation  to  yield  the  number  of  marbles  falling  within  the  rink. 

In  a second  approach,  knowing  the  boy's  distance  from  the  rink,  wire- 
sample  paths  of  marbles  could  be  examined  to  see  if  they  fall  within 
the  rink,  and  to  see  how  long  the  allowed  stopping  distance  is  tor 
each.  This  approach  may  seem  less  direct  compared  to  the  former;  but 
its  results  predict  the  marbles  landing  in  the  ri.nk  for  each  path 
individually,  without  dependence  on  the  width  of  a cup’s  mouth. 

In  the  analogy,  a very  complicated  rink  represents  the  wide-angle 
stopping -muon  telescope?  the  ireirbles  are  rruons;  and  the  cups  model  the 
zenith-angle  detectors.  The  correlation  of  stopping-muon  flux  to 
angular  response  was  to  be  accomplished  using  the  zenith-angle  detectors 
in  coincidence  with  the  muon  telescope  as  a whole.  The  mobility  of  the; 
zenith-anqle  detectors  would  allow  them  to  covex , over  a period  of 
several  months,  the  entire  viewing  aperture  of  the  muon  telescope. 

The  expected  count  rate  frem  such  an  arrangement  was  to  be  predicted 
by  a numerical  integration  scheme  by  ccrputer.  The  rate  predictions 
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of  this  scheme  depended  upon  the  mathematical  description  of  the 
detector  geometry,  just  as  an  est invite  of  the  number  of  marbles 
each  cup  would  collect  would  depend  upon  its  position  within  the  rink. 

The  rate  predictions  for  the  zenith-angle  detectors  were  to  be 
obtained  for  several  flux  distributions.  By  curve-fitting,  the  flux 
distribution  yielding  the  closest  natch  of  predicted  and  observed 
count  rates  was  deemed  to  be  the  extant  distribution  within  the 
surrounding  building.  An  analysis  of  the  wide-angle  telescope's 
response  could  then  be  performed  either  by  numerical  integration  or 
by  the  "trajectory"  method.  The  latter  method  is  actually  a pseudo- 
Monte  Carlo  technique,  vhich  determines  for  nany  particle  trajectories 
whether  or  not  each  meets  the  geometr  ical  criteria  for  detection. 

All  of  these  predictive  processes  involved  iterative  schemes, 
vhich  v*?re  programmed  on  the  Honeywell  635,/Dartmouth  Time-Sharinq 
System  at  the  United  States  Naval  Academy,  in  programs  using  the 
BASIC  language. 

B.  Detector  Response  by  IMnerical  Integration 

Figure  10  depicts  the  two-detector  problem  variables.  For  two 
elemental  detector  areas  dA,  and  dA2 , which  are  separated  by  the 
vector  r12,  and  whose  normals  Nj  and  N2  make  angles  and  02  with 
r j2>  the  absolute  count  rate  seen  by  these  detectors  in  coincidence 
will  be 


dR  =-  I f.l»2  (dA2  cos  *>2)  c5Al  008  °l 


r12 


(1) 
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where  Cj  and  f 2 are  the  detection  efficiencies  of  dAj  and  dA2,  which 
shall  be  assured  to  be  unity  for  the  present.  I is  the  intensity  of 
penetrating  radiation  in  the  direction  rj2.  I norrrally  has  the  units 

particles 

cr»2.Sr.sec 


Ihe  term  dAj  oos  Oi  jUst  tjle  901^  angle  d ft  determined  by  dAj  and 


r12*  Thus,  the  exact  count  rate  of  a two-detector  system  is  given  by 


(r)dA  >pQ 


(2) 


where  I (r)  is  the  intensity  function  of  incoming  particles,  A is  the 
area  of  one  detector,  and  ft  is  the  solid  angle  the  second  detector  subtends. 

The  counted  rate  of  stopping  particles,  however,  is  influenced  by  the 
energy  spectrum  of  the  incoming  particles  as  well  as  by  the  flux  inten- 
sity. A given  particle  mist  remove  all  of  its  momentum  via  seme  process 
(in  the  case  of  a muon,  by  ionization)  in  order  to  stop,  and  thus  must 
transit  a given  thickness  of  material.  For  an  elemental  detector  area 
dAj  and  an  elemental  stopping  thickness  dA^Z,  the  stopping  particle 
problem  is  diagraitmed  in  Figure  11.  Assuming  100%  rejection  of  through 
particles,  the  stopping-particle  count  rate  is 


dR  - I (dA2  cos  <*2)  ***1  ” 91  dp  , 


'122 


(3) 
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where  dt  - 


dZ 

coeTX 


, and  o is  the  stopping  power  of  the  detector 


material,  expressed  as  momentum  loss  per  unit  length.  The  intensity 
1 is  now  the  differential  flux,  expressed  as 


particles 


an  .sr. sec. momentum  unit 


For  the  particular  situation  of  counting  stopping  muons  at  sea  level, 
the  intensity  I is  a slowly-varying  function  of  energy,  and  a value 
of  I is  selected  which  corresponds  to  the  mean  value  of  energy  deposited 
by  a stopping  particle. 

The  oount  rate  of  a stopping-muon  detector  mav  thus  be  expressed  as 


R -o 


E)dA  di-'  dp  , 


(4) 


or,  knowing  the  Cartesian  form  of  the  differentials,  as 


■ -m 


I(r,E)dx2dy2cos  02 


dxidyicos  dz2 


r12‘ 


COS  02 


(5) 


The  tern*  ooe  02  cancel,  and  uie  term  r^2  nay  be  found  if  tne  area  d/.| 
and  the  volume  dA^Z  are  given  Cartesian  coordinates: 
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|r12!  = y (x2-Xl)2  + (y 2-y^' 


(z2 -zi)‘ 


(6) 


Additionally,  cos  0^  may  be  found  if  the  orientation  of  ilfy  is  kncwn 
with  respect  to  dA2.  In  this  experiment,  aAj  is  parallel  to  dA2  for 


the  wide-angle  muon  telescope.  (The  x-y  plane  is  chosen  pared  lei 
to  the  plates  and  the  stop  tank's  top  face) , making 


cos  8}  = 


21  ~ z2 
*12 


For  the  zenith-ancle  detectors,  conversely,  due  to  their  radial  posi- 
tioning,  dAj  is  nearly  perpendicular  to  ri2  for  all  positions,  hencr 


cos  6]_  = 1. 

The  integration  may  new  be  performed  as  soon  as  the  expression  for 
I(r)  is  specified. 

Traditionally,  the  integral  oosmic-ray  distribution  has  been 
expressed  as  a power  of  the  cosine  of  the  zenith  angle,  that  is,  as 

I (©,♦)  - I0  cosn  0 , 

where  IQ  is  the  intensity  of  particles  arriving  vertically,  and  O and 

( 1 , 9 , 9 \ 

♦ are  the  zenith  and  azimuth  angles  of  the  particle's  path.  An 

exponent  of  zero  corresponds  to  an  isotropic  flux,  as  is  essentially 

• • l \ * ) 

the  case  with  the  primary  radiation. ' At  sea  level,  an  exponent 
greater  than  zero  is  generated  fcy  attenuation  throuah  the  greater 
atmospheric  thickness  traversed  along  larger  zenith  ancles,  and  in  the 
case  of  muons,  by  the  greater  decay  probability  along  the  longer  paths 
of  larger  zenith  angles. 

A differential  flux  distribution  I(0,+,E)  may  be  expressed  as  a 
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coerne-pcwer  law  for  any  specific  energy  or  energy’  band.  'Ihe  exponent 
n may  vary  frcr  energy  bo  energy  within  the  measured  range;  however, 
as  stated  previously,  the  slowly-varyinn  magnitude  of  the  flux  with 
energy,  in  the  particles  involved  in  this  experiment,  allowed  tlie 
substitution  of  a sinole  flux  expression  where  the  vertical  flux  value 
was  averaged  ever  the  range  of  energies  involved. 

The  approach  of  the  cosine  power  law  was  adopted  in  this  experiment. 
The  deployment  of  this  function  has  not  always  yielded  a fit  of  high 
accuracy,  but  it  does  reflect  its  usage  in  the  literature. 

Bipioyingi  the  mean  cosine  power  expression  for  the  flux, 
equation  (5)  becomes 


R « <rl 


<Zl-Z2> 


n+1 

• dxjdy^dxjdy-jdZj 


(9) 


far  the  case  of  parallel  detectors,  and  hence  for  the  wide-angle  muon 
telescope.  In  the  case  of  a zenith-angle  detector,  for  vhich  cos  ‘^-1 , 
equation  (9)  becomes 


TJ 


"///// 


(Zi  — Z-i) 

1 * dx^dy^dx^dy^dZ-,  . 


(10' 


Appendix  B contains  a listing  and  a sasple  run  of  each  of  1 HASIC- 
language  programs  designed  to  predict  the  zenith-angle  detector  coin- 
cidence rates  using  equation  (10).  For  consistency,  a grid  spacing  of 
2"  x 2"  x 2"  was  used  in  all  programs. 

Program  ZADRESP  is  predicated  on  the  assumption  that  each  zenith- 
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angle  detector  (hereafter  called  a ZAD)  is  sufficiently  small 
that  it  may  be  treated  as  an  elemental  area;  in  terms  of  equation 
(5) , dA|  is  equal  to  the  entire  SAD  scintillator  area.  For  pro- 
grams 2ADRESP2  and  ZADRESP3#  the  ZAD  is  divided  into  sixteen  smaller 
dAj  s.  For  the  latter  two  programs#  the  coordinates  of  each  dAi  was 
first  written  into  a file  by  the  program  ZADWRTIE#  from  which  file 
it  was  extracted  in  the  course  of  the  numerical  integration. 

The  radial  distance  from  the  stop  tank  center  to  the  ZADs  being 
fixed  at  63.0",  the  position  of  each  ZAD  may  be  specified  by  r,  and  £, 
the  spherical  coordinates  of  its  center  relative  to  the  stop  tank 
center.  These  coordinates  are  different  from  6 and  +#  the  parameters 
of  particle  arrival  path#  since  each  ZAD  counts  particles  with  a range 
of  directions.  Figure  12  shows  the  variety  of  particle  paths  which  a 
single  ZAD  counts,  in  coincidence  with  the  wide-angle  telescope. 

Predictions  made  by  program  ZADHESP3  for  various  flux  descriptions , 
i.e. , for  various  values  of  the  cosine-power  n,  are  shewn  in  Figure  13. 
The  chi-square  technique  to  be  described  in  Section  V ccnpared  these 
predictions  to  the  actual  data  obtained  from  the  ZADs.  Assuming  that 
an  acceptable  confidence  level  is  obtained  from  the  chi-square 
analysis#  it  is  then  valid  to  use  the  detector  description  of  the 
ZADRESP  programs  in  another  program  which  predicts  the  response  of 
the  wide-angle  telescope  as  a whole.  Program  MJRESP,  in  Appendix  C, 
predicts  the  angular  response  for  the  wide-angle  telescope,  as  well 
as  the  half-angle  in  the  north-south  and  east-west  planes.  Its  out- 
put is  in  units  of  Iq  a#  where  a is  the  scintillator  stepping  pewer. 
These  results  allow  the  oenparison  of  the  measured  rate  with  the 
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vertical  flux  measurements  published  by  other  authors,  giving  a 
value  for  detector  efficiency  of  the  stopping-muon  system. 

MLJREKP  will,  in  addition,  calculate  the  count  rate  expected 
within  a given  zenith-angle  aperture.  In  this  manner,  the  aperture 
within  vhich  half  of  the  detected  particles  fall,  the  so-called  hall- 
angle,  may  be  found  directly.  In  the  north-south  plane,  the  halr- 
angle  indicates  the  expected  influence  of  cosmic-ray  sources  to  the 
north  and  south  of  the  zenith  circle  on  the  celestial  sphere.  In 
the  east-west  plane,  the  half-angle  indicates* the  expected  time-width 
of  an  anisotropy  in  cosmic-ray  rate,  in  either  sidereal  or  solar  time. 
Figure  14  shows  the  count  rate  as  a function  of  zenith-angle  aperture 
for  various  flux  distributions.  The  half -angle,  it  may  be  seen,  is 
a weakly-varying  function  of  the  flux  cosine-power. 

C.  Detector  Response  by  the  Trajectory  Method 

Assume  that  a muon  stopped  at  seme  point  within  the  stop  tank, 
and  arrived  along  some  path  specified  by  0 and  $.  One  can  determine 
whether  or  not  the  particle  would  have  been  counted  according  to 
whether  the  path  passed  inside  or  outside  the  boundaries  of  the  top 
coincidence  plate.  Performing  this  operation  for  every  point  within 
the  stop  tank  (or  for  the  uniform  matrix  of  points  used  in  the  ZADRESP 
programs)  will  yield  the  fraction  of  particles  of  each  arrival  direction 
(6,4>)  that  are  counted.  If  the  results  are  sunned  in  azimuth,  the 
summation  yields  the  relative  ability  of  the  nuon  telescope  to  count 
particles  arriving  from  various  zenith  angles.  This  leiativo  ability 
is,  of  course,  independent  of  the  1 l>ix  distribution.  pr<  <:c-:s  \ as 
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performed  only  for  the  wide-angle  telescope  and  not  for  the  zenith- 
angle  detectors,  due  to  the  small  solid  angle  that  each  2AD  subtends. 

Program  DETECTOR,  listed  in  Appendix  C,  performs  this  process, 
ass uning  a symmetrical  detector.  Program  DETBCTR2  uses  more  accurate 
detector  dimensions,  but  yields  the  same  zenith-angle  response  data. 
Figure  15  plots  the  relative  counting  rate  of  the  wide-angle  telescope 
as  a function  of  particle  zenith  angle  of  arrival. 

Multiplying  the  relative  counting  ability  by  the  flux  distri- 
bution for  each  angle,  then  multiplying  again  by  the  zenith-angle  sine 
curve,  yields  the  same  aperture  curve  that  program  MJRESP  does.  The 
DETECTOR  programs  are  valuable  because  they  isolate  effects  of  the 
detector  geometry  frcm  effects  of  the  flux  distribution. 


Detection  ability  relative  t 
for  a vertical  particle 


rigure  15.  Relative  detection  ability  of  wide-angle  tele 
scope,  generated  by  program  DETEJTR2. 
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IV.  BCPSKMWEAL  DCTHWDBffKW  OF  ZENnH-ANQJS  DISTRIBUTION 
A.  OdI  lection  of  Data 

Initially,  a wide-angle  rruon  rate  of  800  counts  per  hour,  and 
a ZAD  solid  angle  of  0.5%  of  the  total  telescope  solid  angle,  were 
assumed.  A statistical  accuracy  of  5%  was  desired,  requiring  400 
counts,  which  should  then  have  been  available  in  100  hours.  To  allow 
for  variation,  counts  were  made  for  one  week  in  each  ZAD  configuration, 
later,  it  was  observed  that  the  100-hour  figure  was  an  under  estimate. 
The  statistical  error  values  will  be  stated  along  with  the  measurements. 

For  each  count,  due  to  limitations  in  the  umber  of  coincidence 
and  discriminator  units  available,  only  three  ZADs  could  be  operated. 

Of  tiie  six  ZAD6  available,  three  were  chosen  far  their  uniformity  of 
response,  and  were  not  replaced  during  the  course  of  the  experiment. 
These  ZADs  were  Timbered  4,  9,  and  15.  ZAD  timber  4 was  optically 
recoupled  once  during  the  course  of  the  experiment,  but  neither  its 
gain  nor  its  response  were  affected. 

Taking  of  data  uu— need  on  22  October  1974,  and  was  acapdeted 
on  28  April  1975,  with  data  being  retrieved  approximately  weekly  fran 
the  automatic  data  output  system.  No  data  was  taken  during  the 
periods  17-26  December  and  14-25  January.  All  data  was  in  the  form 
of  typed  and  paper  tape  teletype  outputs. 

Normalized  count  rates  were  used  rather  than  absolute  rates  in 
the  data  analysis.  Normalization  was  done  relative  to  the  ZAD  rate 
in  the  vertical  position;  thus,  for  each  weekly  count,  one  of  the 
three  ZADs  was  required  to  be  in  the  vertical  slot*  Normalized 


races  were  used  far  the  following  reasons: 

— ooamic  rays,  ^specially  low-energy  ones,  are  susaeptible 
to  varying  attenuation  with  atrtaspheric'  temperature  and  pressure 
variation;  assuaina  that  these  variations  would  affect  the  flux 
equally  at  all  angles,  the  normalization  could  eliminate  pressure- 
and  temperature- linked  variations; 

— the  flux  description  I * 1q  oosn0  is  suited  to  a normalized 
description,  since  normalization  eliminates  dependence  on  the  exact 
value  of  Ig;  and 

— gain  variations  between  tubes  could  be  compensated  for  in 
normalizing,  but  not  for  absolute  count  rates. 

Section  B describes  the  compensation  for  ZAD  efficiency,  and 
Section  C describes  the  calculation  of  statistical  error.  TTva 
sections  following  Section  C explain  the  experimentally-determined 
corrections  for  various  errors. 

B.  Detection  Efficiency  of  Zenith- Angle  Detectors 

It  was  initially  assured  that  each  ZAD  was  so  constructed  that 
its  efficiency  was,  if  not  100%,  very  nearly  equal  to  that  of  each 
other  ZAD.  Preliminary’  examination  of  the  data,  however,  showed 
the  inaccuracy  of  this  assimpticn.  The  method  of  reciprocity  was 
applied  to  compensate  for  the  differing  efficiencies. 

Consider  a zenith-angle  detector  A,  which  has  the  position 
(C,  £) , and  a second  detector  B,  in  the  vertical  position.  The 
rate  counted  by  each  detector  is  equal  to  toe  true  rate  (that 


determined  by  detector  geometry  alone)  multiplied  by  the 


efficiencies  of  each  detector  component;  that  is, 

R^r>  #0 


and 


Rg(0)  = CgE^EgRQ 

where  R(:,)  and  Rq  are  the  true  nuon  rates  through  the  angled  and 
vertical  positions,  respectively.  Assuming  that  the  top  plate  and 
stop  tank  efficiencies  are  independent  of  2enith  angle  (see  Section 
r / F) , the  apparent  normalized  rate  becomes 

RA(t,C)  Ca  R(C,0 

Rfi(O)  ~ Eg  Rq 

Suppose  that  these  detectors  are  new  reversed,  and  a second 
count  is  taken,  with  B at  (r,,Z) , and  A at  the  vertical.  This  apparent 
normalized  rate  is  then 

R ^ Eg  R(C»£) 

R^fO)  ea  Rq 

Taking  the  product  of  the  apparent  normalized  rates,  the  efficiencies 
cancel,  and 

R(C,0  / RAk,C> 
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Thus,  the  efficiency  difference  may  be  otmcmneated  for  without 

Knowing  the  individual  efficiencies,  by  taking  two  counts  with  two 

detectors  in  reversed  locations.  The  quantity  & will  be 

RO 

referred  to  as  a reciprocity-nornelized  count  rate. 

Table  4 lists  the  reciprocity-normalized  rates  for  all  of  the 
mount  positions  at  which  measurements  were  taken.  The  second  column 
contains  the  raw  mater  of  counts  in  each  detector,  the  vertical 
count  being  the  bottom  number  in  each  case.  Each  ratio  was  con- 
ducted with  the  same  set  of  tubes  as  its  reciprocal,  that  is,  tubes 
4 and  9,  9 and  15,  or  15  and  4,  with  the  exception  of  the  data  at 
(t,£)  « 20*  north,  where  all  three  tubes  were  used,  and  reciprocity 
was  obtained  at  that  point  with  the  relation 


R(20#N)  Rg(20)  %5(20)  1^(20) 


The  third  col  ism  of  the  table  contains  the  reciprocity-normalized  rates. 


C.  Statistical  Error  Calculation 


The  arrival  times  of  cosmic  ray  particles  may  generally  be  treated 
as  random  occurrences.  If  a rate  is  indeed  generated  fy  randomly- 
distributed  counts,  the  standard  deviation  of  a nutter  of  counts  ray 
be  conservatively  estimated  as  the  square  root  of  that  muter;  that  is, 


I ! 


A measured  normalized  rate  may  be  expressed  as 


A t A 


R(;,u  _ 

R(0)  B i & 


where  A and  B counts  are  obtained  by  two  detectors  in  the  same  time 
period.  If  this  fraction  is  expanded,  the  deviation  generated  is 
probably  somewhat  larger  than  the  standard  deviation,  a result  of  the 
propagation  of  errors.  The  expansion  is 


-1 


A ■■  = (A  i v'ft-d  i if) 

B i B B 


i An  + /A  , , 

gU  - -g  i -g)  • 


where  second-order  terms  have  been  dropped.  A conservative  statement 
of  the  standard  deviation  of  a normalized  rate  is  then 


r * , 1 + J.  . 

» bI  ' 


For  a reciprocity-normalized  rate,  given  that  A and  C are  the 
numbers  of  particles  counted  at  a particular  mount  position,  and  B 
and  0 are  the  vertical  counts  for  the  corresponding  time  interval. 
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(18) 


(19) 


(20) 


J 
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the  randcaness  error  gives 


(21) 


again  neglecting  second-order  teems.  An  exaggerated  estimate  of  one 
standard  deviation  is  given  when  all  error  terms  in  equation  (21)  are 
either  positive  or  negative,  i.e. , 


This  is  actually  the  unbiased  estimate  of  crp^,  rather  than  the  ? 
generated  by  data  distribution.  The  last  oolixan  of  Table  4 contains 
the  Cp»4  for  each  data  point* 

All  mathematical  eliminations  and  error  predictions  have  now  been 
performed.  Experimental  determination  of  certain  systematic  errors 
is  explained  in  the  following  sections. 


n.  Shower  Corrections 

Often,  following  a primary  particle's  interaction  with  an  atmo- 
spheric nucleus,  the  reaction  products  have  sufficient  energy  to 
cause  further  nuclear  interactions.  The  propagation  of  such  a 
shower  of  saoondarv  cosmic  rays  is  such  that  nmnv  particles  may 


4 


arrive  at  the  earth's  surface  simultaneously.  Should  a stopping 
nuon  be  associated  with  other  particles  in  a shower,  a zenith-angle 
detector  may  be  erroneously  triggered  by  an  associated  particle, 
the  resulting  event  signature  mimicking  that  of  a muon  passing 
through  the  ZAD.  Figure  16  demonstrates  the  generation  of  a shower- 
produced  count. 

Because  the  arrival  directions  of  showers  are  strongly  concen- 
trated about  the  vertical,  ZADs  at  mount  positions  nearer  the 
vertical  were  not  expected  to  accrue  shower-generated  aounts,  since 
many  of  the  associated  particles  should  pass  through  the  bottom  ant  i- 
od incidence  plate,  negating  the  entire  event.  Detectors  with  ’ = 40° 
in  any  plane,  and  with  * * 30°  in  the  north-south  and  east-west 
planes,  however,  did  not  project  onto  the  bottom  plate,  and  were 
thus  deemed  to  be  susceptible  to  shower-generated  counts. 

A means  of  correcting  the  data  from  these  points  was  created  by 
mounting  a detector  vertically  beneath  the  scintillator  of  a 10"  or 
40°  detector.  In  this  position,  tfuwn  in  figure  17,  the  lower 
detector  is  susceptible  to  approximately  the  same  number  of  shower- 
generated counts  as  the  7AD  above  it;  however,  because  it  is  mounted 
below  the  plane  of  the  top  coincidence  plate,  it  can  register  no  true 
counts.  The  reciprocity-normalized  rates  of  these  shower  detectors 
were  applied  as  negative  corrections  bo  the  corresponding  ZAD  rates. 
The  shower  correction  was  found  to  be  (-.024  * .006) Rq  at  \ * 40°  , 
and  (-.064  t .01)Rq  at  z,  ■ 30*. 
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Figure  1f>.  Shower-Generated 
ZAD  count  (not  to  ec^le.) 


Figure  17.  ucc-tion  : f 
Jetector  for  exp^ri^*  :.t  » . 
ietermirr  tion  of  n.ow<  r 
correction,  f « ^cf  (no* 
to  eciIp,! 


- — -V»  , 


*^**-*“  - '***  • 


E.  Gorrelated-Pulsing  Corrections 


Due  to  the  geanetry  of  the  wide-angle  muon  telescope,  it  is 
possible  for  sane  through  particle  trajectories  to  trigger  a TSB 
count  by*  missing  the  bottom  plate,  as  shown  in  Figure  18.  Further, 
some  of  these  trajectories  may*  intersect  the  boundaries  of  a ZAD  in 
the  (f.,0  * 40°  north,  east,  south,  and  west  positions.  Should  a nv 
S pulse  be  generated  in  the  acceptance  interval  thereafter,  that  is, 
between  0.3  and  5.3  nsec  later,  an  erroneous  count  will  be  generated. 
Hie  SfE  pulses  following  a false  TSB  count  nay  acme  either  frars  the 
random  inpact  of  aosmic  rays,  from  randan  noise,  or  from  a correlated 
second  pulse  generated  within  the  system. 

A multichannel  analyzes:  accumulation  of  the  "stopping-muon" 
decay  times  was  available  for  all  data  taken  between  22  October  1974, 
the  beginning  of  this  experiment,  and  1 February*  1975.  Gonparison 
of  this  time-dona  in  data  with  the  theoretical  muon  decay  curve  indi- 
cated that  approximately  20%  of  all  ttuon  counts  accumulated  to  that 
time  were  generated  by  either  correlated  pulsing  or  accidental 
coincidences,  lb  test  the  possibility  that  accidental  coincidences 
were  responsible,  the  TS§  and  ? T B rates  were  measured,  and  the 
accidental  rate  was  calculated.  The  results  were: 


counts 

R-ISB  - 175  -SGT  * 2.92  cps 
counts 

Rsfg  - 2100  "nUJi — * 35  cps 


Tr  * decay  pulse  acceptance  window 
« 5.0  x lO1"6  sec 


no  B p'ueo 

produced 


Figure  18.  Generation  of  TSB  pulses  by  through  particles, 
leading  to  acceptance  of  false  auon  signatures. 
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These  figures  Yield  an  accidental  rate  of  3.11  x 10“*  ops.  The 

raw  muon  rate  at  that  time  was  anproximately  900  gP"-ints,.  or  . '*3 

hour 

cps,  indicating  that  the  20%  sxcess  in  counts  oould  not  have  been 
generated  by  purely  randan  coincidences. 

At  least  one  known  phenomenon  can  explain  second  pulses  in  the 
stop  tank  that  are  coupled  to  those  produced  by  cosmic  rays.  Should 
a gas,  such  as  air,  leak  into  one  or  both  of  the  5"  photomultiplier 
tubes,  its  ionization  by  the  accelerated  photoelectrons  from  a 
scintillator  flash  would  make  it  a charge  carrier.  Sufficient  quan- 
tities of  such  charged  particles  dislooing  note  electrons  from  the 
photo-cathode  could  cause  one  or  many  correlated  pulses. 

Correlated  pulsing  was  observed  on  an  oecillosaope  when  the 
5“  tubes  were  decoupled  from  the  stop  tank  and  triggered  bv  a pulsed 
light  source.  The  correlated  pulsing  was  at  all  times  less  than 
130  millivolts  in  amplitude,  and  was  concentrated  about  a 0.8  usee 
lag  behind  the  initiating  pulse.  This  lag  corresponded  to  a pro- 
minent peak  in  the  time-domain  spectrin. 

Approximately  80%  of  the  correlated  pulsing  was  eliminated 
upon  reducing  the  5"  phototube  voltage  from  2300v  to  2200v.  After 
recoup ling  the  5"  tubes  to  the  stop  tank,  the  correlated  pulsing 
was  further  reduced  by  increasing  the  stop  tank  discriminator  thresh- 
olds from  20  mv  to  100  mv.  A variety  of  discriminator  thresholds 
was  examined,  100  mv  being  the  lowest  possible  setting  which 
coincidences  between  through  aosmic  rays  (TBS)  and  stop  tank  second 
pulses  (STB). 
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Following  this  correction,  the  average  stopping  muon  count 
rate  decreased  from  about  900  counts  per  hour  to  about  540  counts 
per  hour.  This  40%  decrease  was  somewhat  greater  than  the  drop 
expected  from  the  elimination  of  aorrelated-pulse  counts,  implying 
that  the  detection  efficiency  of  the  system  had  decreased.  However, 
the  tjme-danain  data  now  showed  no  noise-generated  peaks  or  other 
anomalies,  and  the  half-life  obtained  from  it  was  within  2%  of  the 
published  nuon  half-life  in  organic  compounds.  * 

However,  a time-domain  accumulation  of  "muons"  detected  by  the 
40°  S detector  only  showed  that  approximately  20%  of  these  counts  were 
still  caused  by  multiple  pulsing.  This  is  not  inconsistent  with  a 
oorrelated-pulsing  rate  of  0.3%  or  less  from  the  wide-angle  telescope 
as  a whole.  The  requirement  was  set  that  no  data  from  t)«  40°  N or 
40°  S positions  be  accepted  from  the  time  period  before  the  reduction 
of  the  stop  tank  tube  gain,  lb  the  accepted  data  was  applied  a 
correction  of  (-.20  ± *05)R4on,S,  or  {-.020  ± .005)Rq. 

F.  Compensation  for  Counting  Efficiency  of  Muon  Telescope  Opponent* 

As  has  been  shewn  in  Section  IV  A,  the  detection  efficiency  of 
the  zenith  angle  detectors  can  be  corrected  for  bv  switching  two 
detectors  between  an  angled  position  and  the  vertical.  However,  nr 
such  sinple  means  exists  for  correcting  for  variations  in  efficiency 
within  the  large  acmponents  of  the  muon  telescope,  al thcuqh  these 
efficiency  variations  are  measureable. 

The  efficiencies  at  various  positions  on  the  top  and  bottom  plates 
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were  measured  by  the  "sandwich"  method,  using  the  detectors 
described  in  Section  II  C.  Figure  19  is  a schematic  representation 
of  two  snail  detectors  being  used  to  measure  a local  efficiency  of 
a larger  planar  detector  P.  Sn^  and  Sr^  represent  the  smaller 
"sandwich"  detectors.  The  count  rate  of  Sn^  and  Sn2  in  coincidence  is 


r12  = elc2  GI0 


(21) 


vhere  and  e2  are  the  sandwich  detector  efficiencies,  G is  a 
geometrical  factor  due  to  the  areas  and  arrangement  of  Sn^  and  Sn2» 
and  Iq  is  the  vertical  flux.  It  can  be  seen  frcm  the  figure  that 
every  particle  that  passes  through  both  Sn^  and  Sn2  must  also  piss 
through  P;  the  count  rate  of  Sn^,  P,  and  Sn2  in  coincidence  is  then 
given  by 


r1P2  * clePe2  GI0  ’>  (2'»> 

therefore, 

gp  = . (25) 

R12 

Measurements  of  this  sort  were  conducted  in  various  regions  of  the 
top  and  bottom  plates,  and  the  results  are  recorded  in  Figures  20  and 
21.  These  measured  efficiencies  are  accurate  to  tlO%. 

Figure  20  indicates  that  the  comer  areas  of  the  top  plate  are 
significantly  less  than  100%  efficient.  Particles  detected  by  ZADs 
in  the  40*  NE  and  40°  SE  positions  must  pass  through  these  regions 


Figure  19.  Measurement  of  the  efficiency  of  a planar 
detector  by  means  of  the  "sandwich'’  method.  All  par- 
ticles passing  through  both  and  Snp  are  constrained 
to  pass  through  F as  well. 
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enroute  to  the  stop  tank.  Thus  the  apparent  rates  seen  by  ZADs 
in  these  positions  are  lower  than  the  true  rates.  'The  exact  amount 
of  enhancement  retired  is  uncertain,  Kit  is  probably  between  zero 
and  15».  The  correction  due  to  low  corner  efficiency  was,  then, 

{*-08  * •08)R4ojc#?e>  or  <+‘021  - .021)^. 

It  is  also  apparent,  from  Figure  21,  that  the  lev  efficiency 
regions  of  the  bottom  plate  allow  sane  false  rsft  counts  to  be 
generated.  The  correct ion  for  this  effect  is  inherent  in  the 
correlated  pulsing  corrections  bo  the  40°  N and  40°  S positions  (see 
Section  rr  E) , and  no  additional  correction  is  remired. 

Because  of  the  difficulty  of  measuring  local  efficiencies  within 
the  stop  tank,  and  because  the  snail  stop  tank  volime  was  viewed  b” 
two  high-ef f iciency , high-gain  photomultiplier  tubes,  the  stop  tank 
counting  efficiency  was  assumed  to  be  100%  throughout  the  course  of 
the  experiment. 

G.  Positional  Error  of  Zenith- V«le  betectors 

It  was  felt  that  the  position  of  each  ZAD  was  known  to  t0.5" 
in  any  direction.  This  error  is  insignificant  for  all  except  the 
40°  !1  and  40°  S locations;  as  may  be  seen  in  Figure  13,  a small  chnige 
in  the  reed  position  of  a 40°  N or  40°  S detector  causes  a large 
fractional  change  in  count  rate,  taonrdingly,  after  all  other 
corrections  were  made  to  these  points , 7 4qN  and  ^ 40lq  were  scaled  up 
by  a factor  of  2.S,  due  to  rate  variations  with  position. 
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V.  COMPARISON  OF  MEASURED  DATA  AND  PREDICTS)  RESPONSE 
A.  The  Chi-Square  Fit 


A standard  method  of  determining  which  of  several  theoretical 
expressions  best  fits  a set  of  experimental  data  is  the  "chi-square.! 
test."  In  this  process,  the  quantity 


s-t  V’2 

i-1 

is  to  be  minimized,  where  the  are  observed  data  points,  and  the 
are  the  corresponding  theoretical  values,  selected  from  ary  of  several 
predictions.  In  addition  to  yielding  the  best-fitting  theoretical 
curve,  the  process  can  measure,  although  subjectively,  the  goodness 
of  fit  of  any  theoretical  curve.  The  number  of  points  k,  assuming 
that  the  are  independent  of  one  another,  is  kncr«i  as  the  number  of 
degrees  of  freedom. 

The  magnitude  of  S,  as  expressed  in  equation  (26) , is  applicable 
only  when  the  and  N^  are  total  mntoers  of  counts,  whereas  this 
experiment  deals  with  normalized  rates.  Tb  relate  the  two,  assume 
that  an  observed  rate  R|  is  equal  to  , and  similarly  that  a 

n 

N* 

theoretical  rate  P.  is  equal  to  — . Then 
x Ti 


S * 


i«l 


(27) 


For  a normalized  rate,  R^  trust  be  divided  by  the  normalizing  rate  R^, 

0 

and  by  the  theoretical  rate  P^,  v#iich  cannot  be  expressed  exactly 
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in  terms  of  equation  (27) . An  approximation  may  be  made  by  dividing 
both  the  % and  by  Rj|,  which  gives 


(28) 


The  products  are  simply  the  number  of  aounts  in  the  vertical 
detector  over  the  time  periods  T^.  Equation  (28)  is  also  applicable 
to  reciprocity-normalized  rates,  if  the  substitutions 


(29) 


and 


(30) 


(31) 


are  made,  where  the  A^  and  are  the  raw  numbers  of  ZAO  events  at 
(Cj^,Ci)  > and  the  and  are  the  events  in  the  vertical  ZADs,  for 
each  time  T^.  The  p£  are  simply  the  predictions  of  normalized  rates 
of  a program  such  as  2ADRESP3.  S is  now  used  in  this  form, 

k 

i=l 
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in  the  program  RESULTS,  listed  in  Appendix  D.  This  value  of  S is 
labeled  "FLAT  S"  in  the  output  of  RESULTS,  far  reasons  which  will 
be  apparent  later. 

The  drawbacks  of  the  forma  1a  for  S developed  in  equations  (26) 

through  (32)  are  that  it  is  an  approximation,  improperly  normalized 

as  of  equation  (28) , and  gives  equal  weight  to  each  data  point. 

2 

Equation  (26)  is  actually  valid  only  when  , i.e.,  when  the 

Qi  cure  Russian-distributed  about  the  N£.  Should  sane  of  the  measured 
rates  have  large  known  systematic  errors,  S would  become  unreasonably 
large. 

A more  universally  applicable  expression  for  S is  given  by 


S 


(33) 


where  the  a ± are  the  standard  deviations  corresponding  to  the  Qj 
(after  Wolberg. 1 ? ) The  reasons  for  preferring  equation  (33)  may 
be  stated  briefly.  Using  this  expression,  the  statistically  poorer 
points  receive  less  wei^it  in  the  calculation  of  S.  Further,  any 
known  errors  which  are  systematic  in  nature,  but  vhich  cannot  be  made 
an  exact  correction  to  the  data  (e.g. , the  low  efficiency  of  the  top 
plate  oorner  regions) , may  be  stated  as  an  increase  in  the  value  of 
the  uncertainty  0 ^ . Finallv,  equation  (33)  is  not  limited  in  its 
applicability  to  nurfcers  of  counts;  rather,  the  0^  and  may  be  used 
directly  as  normalized  rates  and  reciprocity — normalized  rates.  Program 
RESULTS  outputs  the  value  of  equation  (33)  as  "WEIGHTED  S". 


i 
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The  actual  use  of  the  variable  is  as  a reference;  that  is, 
values  of  x are  tabulated,  as  a function  of  certain  "confidence 
levels,"  the  measure  of  goodness  of  fit  mentioned  earlier,  and  of 
k,  the  number  of  degrees  of  freedom.  S is  compared  to  these  tabu- 
lated values. 

B.  Analysis  of  Data 

Table  5 is  a suaoary  of  all  of  the  corrections  in  Section  TV, 
along  with  the  net  results  upon  each  normalized  rate  and  e . These 
were  the  corrections  employed  in  program  RESULTS,  which  performed  the 
S-minimi ration  process  by  evaluating  S for  various  values  of  the 

a 

cosine-power  n.  In  Figure  22,  the  weighted  S values  generated  by 
equation  (33)  are  given  for  various  values  of  the  cosine-power.  The 
upper  curve  was  generated  using  uncorrected  values  for  the  reciprocity- 
normalized  data  points,  and  is  shown  for  comparison  only.  The  lower 
curve,  when  all  corrections  have  been  applied,  clearly  shows  that  S 
is  minimized  when  the  cosine-power  is  about  4.7.  Further,  published 
tables* 1 7 'P* * ^ nay  be  used  to  find  the  relative  probability  that 
the  data  came  from  each  set  of  predicted  values.  Frcm  these  relative 
probabilities,  the  "confidence  levels"  described  previously,  it  nay 
be  seen  that  n probably  deviates  frcm  the  value  4.7  by  no  more  than 


10.5: 


n 

weighted  S 

3.5 

25.5 

4 

13.1 

4.5 

7.64 

5 

7.77 

5.5 

12.0 

Relative  Closeness  of  Fit 

<0.01 

0.2 

0.7 

0.7 

0.3 


• •’SWSJs’S 
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In  order  to  show  the  influence  of  the  marginal  data  points 
on  the  value  of  S,  the  curves  in  Figure  23  were  generated  by  data 
excluding  that  at  the  marginal  points  40°  N and  40°  S.  Cdgaro  the 
clotted  curve  to  the  uncorrect eri  1 Q-poi nt  curve  at  the  top  of  the 
figure.  The  dotted  curve,  which  excludes  the  marginal  points  taut 
employs  no  other  corrections , shows  that  the  error  m the  uncorrected 
data  is  due  mainly  to  the  40°  N and  40°  S data.  Adding  the  corrections 
to  the  8-point  curve  generates  the  solid  curve  which,  it  may  be  seen, 
is  very  similar  to  the  10-pomt  curve  in  Figure  22. 

The  flat-S  values  generated  by  equation  (32)  are  plotted  in 
Figure  24.  They  were  not  employed  in  the  data  analysis,  because 
equation  (32)  is  not  considered  to  weight  the  data  points  proper ly . 

The  confidence  levels  quoted  in  this  section  should  not  be  inter- 
preted as  a measure  of  confidence  of  the  observer.  Rather,  the> 
represent  a relative  measure  of  the  closeness  of  fit  of  each  set  of 
predicted  to  the  measured  data. 
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Figure  25.  Comparison  of  S for  data  exclusive  of  i*(f  !*  and 
W 3 data  (k  » 8 degrees  of  freedom),  with  S curve  gen- 
erated by  10  uncorrected  data  points  (k  » 10).  Computa- 
tions by  program  RESULTS. 
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Figure  2 4.  Flat  (unweighted)  L curvet,  for  all  data  (k  - 
10  degrees  of  freedom),  generated  by  uncorrected  and 
corrected  data.  Computations  by  program  RESULTS. 
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VI . DISCUSSION  OF  RESULTS 
A.  Mathematical  Description  of  Detector  System 

Ihe  program  MURESP  is  an  exanple  of  inductive  logic  --it  springs 
from  the  description  of  the  zenith-angle  detection  system  in  the 
ZADRESP  programs.  Its  validity,  therefore,  depends  upon  two  cri- 
teria of  induction:  the  accuracy  with  which  the  ZADRESP  programr- 

predicted  the  response  of  the  wide-angle  telescope  in  coincidence 
with  the  zenith-angle  detectors;  and  the  validity  of  replacing  the 
zenith-angle  detectors  with  the  top  coincidence  plate  as  a whole. 

Figure  25  shows  the  corrected  data  in  the  north  lialf-plane  and 
at  the  40°  northeast  position,  superimposed  upon  the  ZADRESP3-generated 
rate  predictions  for  n » 4.5  and  n = 5.  The  error  bars  are  of  the 
values  ocaputed  by  program  RESULTS.  All  corrected  data  points  are 
within  one  standard  deviation  of  the  predicted  curve,  indicating  that 
ZADRESP3  indeed  provides  an  accurate  description  of  the  ZAD-ooupled 
wide-angle  telescope  response,  satisfying  the  first  criterion  for 
induction. 

It  is  assumed  that  the  second  criterion  is  satisfied,  since  the 
area  elements  of  the  top  plate  tray  be  made  arbitrarily  small.  The 
resulting  detector  description  is  then  identical  to  that  fear  a ZAD 
except  for  a factor  of  cos  0.  The  predictions  of  MURESP  may  then  be 
said  to  be  accurate  within  the  limits  of  this  experiment.  These  pre- 
dictions are: 

— a wide-angle  telescope  half-angle  (zenith-angle  within  which 
half  of  all  particles  are  detected)  of  23.0°  ♦ .25°,  or  a full  width 


n • 5 


n * i+.5 


Cf  Kf  2(J  5#  -»tf 

Detector  zenith  angle,  0 

Figure  25 • Comparison  of  predictions  fro-'  promr 
to  measured  rotes  in  north  half-plane  (curve  e*tr«p 
ted.  to  4Cf  NE),  for  n -*  4.5  and  n » 5.  Cslculatio. 
rates  from  raw  data  by  program  RLEUi.T.',. 
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at  half-maximum  of  46.0°  + 0,5°,  taking  the  value  of  n to  be 
4.7  + 0.5;  and 

— an  average  zenith-angle  secant  for  the  particles  detected 
of  1.19,  which  is  important  in  energy  computations. 

3.  Description  of  the  Stopping-Muon  Flux  at  Sea  Level 

The  value  of  the  cosine-power  n for  stopping  muons,  as  obtained 
in  this  experiment,  is  significantly  higher  than  the  cosine-powers 
of  the  differential  muon  energy  spectra  observed  by  other  experi- 
menters . For  exanple,  Kraushoar,  in  1949,  measured  the  muon  differ- 
ential range  spectrum  for  R >71  gm  csn"2  of  air  (corresponding  to  a 

\ 

kinetic  energy  of  about  350  MeV) , and  found  n to  be  approximately 
3.3.^'^  Zar,  in  1951,  measured  the  differential  range  spectrum  for 
through  particles  in  approximately  the  same  energy  band,  and  found  n 

( 1 4 ) 

to  be  2.97.  Moroney  and  Parry,  in  1954,  also  obtained  an  n of 
about  3.3  for  the  differential  spectrun  at  0.3  BeV/C. By  contrast, 
the  cosine-pcwer  of  stopping  muons  as  determined  in  this  wor  k was 
4.7  + 0.5. 

Several  factors  contribute  to  this  difference.  First,  all  of 
the  other  observers  cited  above  used  data  as  far  away  from  the  zenith 
as  0 = 60°,  whereas  in  this  experiment  no  measuranents  beyond  *»  = 40° 
were  taken.  The  spectra  in  each  of  the  cited  works  tends  to  dip  until 
the  zenith  angle  exceeds  o = 30°,  when  a leveling-off  occurs;  the 
average  cosine-pcwer  ever  the  entire  range  of  zenith-angles  is  about 
3 to  3.3,  but  the  n in  the  regions  nearer  the  zenith  appears  to  be 


somewhat  greater.  Second,  the  angular  apertures  of  all  of  the 
telescopes  in  the  cited  works  were  somewhat  wider  than  those  of 
the  ZADs  in  this  experiment.  Third,  all  of  the  telescopes  used 
in  the  cited  works  were,  to  sane  degree,  susceptible  to  enhance- 
ment by  vertical  showers  as  the  telescopes  were  inclined  to  larger 
zenith  angles.  The  latter  consideration  nay  be  shown  to  be  insigni- 
ficant, however.  When  shower  corrections  were  applied  to  the  experi- 
mental data  herein,  ti  alue  of  n at  which  S was  minimized  changed 
by  less  than  0.2;  also,  the  cited  experiments  were  consistent  in 
their  values  of  n independent  of  their  shower  susceptibility. 

A fourth  consideration  in  ccnparing  the  zenith-angle  distri- 
butions is  the  energy  band  of  the  particles  bei  lg  counted.  It  is 
known  that  the  cosine-pcwer  increases  with  decreasing  energy,  and 
that  a stopping  muon  is  of  the  minimum  detectable  energy  at  the 
earth's  surface.  The  surrounding  overburden  also  has  an  effect  upon 
the  energy  distribution  at  the  detector.  More  will  be  said  about 
energy  effects  in  the  next  section. 

It  is  reasonable  to  state,  however,  that  the  flux  description 
determined  by  this  experiment, 

1(0)  = IQ  cos4,7*'  , 

is  the  relevant  flux  to  enploy  in  calculations  concerning  the  wide- 
angle  muon  telescope  des'TibU  herein. 

C.  Energy  Considerations 


For  a particle  to  stop  in  the  detection  system,  it  must  pass 


through  three  major  contributors  to  ionization  energy  loss:  the 

building  roof  and  structure;  the  load  soft-ocaponent  shield;  and 
the  scintillator  of  the  stop  tank.  Figure  26  shows  a cross-section 
of  the  building  in  which  the  detection  system  was  maintained.  For 
each  of  the  paths  shown,  the  average  particle  energy  was  calculated, 
using  the  following  assunptions: 

— the  muon  specific  energy  losses  in  scintillating  plastic  and 
in  carbon,  expressed  as  MeV-cm^/gm,  are  approximately  the  same; 

— the  muon  specific  energy  losses  in  concrete  and  in  copper 
are  approximately  the  same;  and 

— no  scattering  is  assumed. 

The  results  of  these  energy  estimates  were: 


Path 

Energy 

Lass 

Zenith  Angle 

Contribution 

\E 

0° 

20  g/cm2 

plastic 

70  MeV 

55  g/am2 

lead 

60  MeV 

90  g/cm2 

concrete 

130  MeV 

TOTAL 

260  MeV 

30° 

23  g/cm2 

plastic 

80  MeV 

63  g/cm2 

lead 

70  MeV 

104  g/cm2 

concrete 

150  MeV 

TOTAL 

300  MeV 

40° 

28  g/cm^ 

plastic- 

100  MeV 

71  g/cm2 

lead 

80  MeV 

450  g/cm2 

concrete 

670  MeV 

TOTAL 

850  MeV 

The  high  energy-loss  values  in  plastic  are  due  to  the  sharp  increase 


Figure  <?6.  Stopring-muor.  telescope  location  within  surround- 
ing building  (CMR  wilding,  Jreenbury  Toint,  U.S. Naval 
Radio  Transmitter  .t.tion,  Annapolis,  Maryland.)  A p:y- 
j/ood-and-shingle  roof  covered  the  area  betAr^n  the  two 
beveled  columns  shown,  which  are  actually  cross-sections 
through  an  octagonal  structure  3.0  feet  thick.  Stopping 
power  of  concrete  blocks  through  which  vertical  cosmic 
rays  must  pass  may  have  been  somewhat  overertiaated , 
since  they  are  porous  blocks,  as  opposed  to  the  dense 
structural  concrete  elsewhere.  Trajectories  of  A Cf  N and 
S actually  transit  different  amounts  of  concrete. 
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in  ^ as  the  muon  stops.  Range  and  energy-loss  values  were  taken 
from  High -Ehergy  Particle  Data,  published  by  Lawrence  Radiation 
Laboratory. 

At  first  glance,  the  nuch  higher  energy  required  to  produce  a 
stopping  muon  with  a zenith  angle  of  40°  seems  to  inply  that  the 
surrounding  building  produces  significant  attenuation  in  the  muon 
rate  at  •>  - 40°.  This  would  tave  artifically  lowered  the  apparent 
cosine-pcwer  n.  However,  upon  examining  the  differential  sea- level 
muon  spectra  of  Moroney  and  Parry  ^ ^ , Kraushaar  * 1 , and  Allkofer 
et  al.  there  is  no  consistent  difference  in  the  differential 
intensities  of  muons  of  energy  260  MeV  and  those  of  860  MeV.  The 
spectra  are  very  flat  in  this  region,  the  greatest  difference  in  the 
two  intensities  being  about  five  percent.  Thus,  the  effect  of  the 
building  as  shielding  is  merely  to  shift  the  energies  of  stopping 
muons  to  another  portion  of  the  energy  spec  trim,  but  not  to  decrease 
their  number  significantly. 

D.  Suggestions  for  Fur tlier  Research 

The  following  possibilities  for  further  research  were  suggested 
by  the  results  of  this  experiment  and  the  problems  associated  therovith. 

First,  it  is  suggested  that  the  zenith-angle  distribution  of 
stopping  nuons  be  measured  at  larger  zenith  angles.  The  only  distri- 
bution of  interest  in  this  experiment  was  that  which  produced  counts 
in  the  wide-angle  muon  telescope.  However,  any  shallowing  of  the 
distribution  at  larger  angles  would  be  of  great  interest  theoretically. 
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Second,  an  expression  for  the  zenith-angle  distribution 
should  be  developed  which  is  a result  of  theoretical  considerations. 

As  stated  previously,  the  use  of  the  expression  I = Iq  oosn*’  is  a 
cus ternary,  rather  than  a theory-based,  procedure.  A prediction  of 
zenith-angle  distribution  should  be  associated  with  each  nuon  pro- 
duction-attenuation model . 

Third,  a more  extensive  study  of  the  incidence  of  si yower -asso- 
ciated stopping  raucns  should  be  undertaken.  Most  showers  in  tins 
apparatus  should  have  been  eliminated  by  the  anticoincidence  plate; 
however,  their  persistence  (due  at  least  in  part  to  the  low  efficiency 
of  the  bottom  plate)  may  permit  a study  to  be  irate  of  them.  It  is 
possible,  also,  that  showers  have  distorted  the  apparent  zenith-angle 
distribution  in  seme  way  other  than  that  assumed  herein. 

Fourth,  a Monte  Carlo  calculation  of  scattering  effects  in  the 
building  and  load  shielding  should  be  undertaken.  Muons  of  the 
energies  being  measured  here  are  at  least  as  susceptible  to  scattering 
as  the  through  particles  measured  in  other  experiments,  for  which 
scattering  corrections  are  often  significant. 

Finally,  two  areas  of  research  have  suggested  themselves  and 
are,  althou^i  not  directly  related  to  zenith-angle  distribution, 
infinitely  related  to  the  meaning  and  importance  of  stopping -muon 
aomnic  rays  at  sea  level.  The  first  is  the  performance  of  a time- 
correlation  of  stopping-muon  arrivals  over  an  acceptance  time  of 
about  one  second.  The  ramifications  of  any  non-randan  time  peaks 
which  might  appear  are  manifold.  Second,  and  lastly,  a collection 
of  shower-associated  stopping-muon  data  should  also  be  undertaken, 


'****—■ 


in  solar  and  sidereal  time.  An  antiquity  exists  in  the  primary 
energy  spectrum  between  lO^eV  and  lO^eV,  and  shower-associated 
stopping  nuons  are  in  the  ideal  range  for  studying  this  ambiguity. 
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Element 

TABLE  1 

Nuclear  Abundances  in  the  Primary 
Comic  Radiation 

Relative  Abundance 

H 

100 

He 

15 

Be 

0 to  0.4 

C,N,0,F 

1.2  + 0.4 

He 

0.2 

*S 

0.09 

Si 

0.07 

Fe,Co,Ni 

0.06 

beyond  Ni 

in 

I 

© 
*— * 

V 

! *!  «Nr*V  *V  *■  ■Jr''l^{’ Vl 


I 
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TAP  E 2 


Major  Dimensions  of  Detector  System 


Size  of  top  and  bottom  plates 

Separation  of  plates 

Scintillator  size  in  stop  tank 

Height  of  stop  tank  center  from 
bottom  plate  upper  surface 

Offset  of  stop  tank  center  from 
centerline  of  plates 


48.0“  x 48.0* 


16  3/4“  x 16  3/4"  x 15“  (N-S) 


15.75" 


0.5"  to  the  south 


Zenith-Angle  Detectors  and  Mount 


Size  of  ZAD  scintillator 


2.0"  x 5.0"  dia. 


Radius  of  mount 


Radius  of  travel  of  ZAD 
scintillator  centers 


63.0" 


Center  of  mount  arc  fran 
center  of  stop  tank 


0 +0.5"  vertical, 

0 + 0.5"  E-W, 

0.?"  + 0.5"  to  the  south 


Height  of  vertical  ZAD  scinti- 
llator frcm  top  plate  upper 
surface 


TABLE  3 


ZMX85P3  Predictions  of  Nomaliaad 
ZAD  Stappinq-Mnon  Count  Rates* 


n,  in  I = Iq  oosn  6 


Z AD 

Position 

2.5 

3.0 

3.5 

4.0 

4.5 

5.0 

5.5 

10°  N 

.9597 

.9519 

.9441 

.9365 

.9290 

.9216 

.9143 

10°  S 

.9597 

.9519 

.9441 

.9365 

.9290 

.9216 

.9143 

20*  N 

.8498 

.8231 

.7974 

.7726 

.7488 

.7258 

.7037 

30*  N 

.6891 

.6407 

.5960 

.5547 

.5165 

.4812 

.4485 

30°  E 

.6857 

.6380 

.5940 

.5533 

.5157 

.4809 

.4488 

30°  S 

.6764 

.6286 

.5844 

.5436 

.5059 

.4710 

.4388 

40°  N 

.0456 

.0379 

.0315 

.0262 

.0218 

.0181 

.0151 

40°  NE 

.5008 

.4384 

.3842 

.3371 

.2961 

.2604 

.2292 

40*  SE 

.4929 

.4313 

.3778 

.3313 

.2908 

.2556 

.2249 

40*  S 

.0232 

.0191 

.0158 

.0131 

.0108 

.0090 

.0074 

* The  program  ZADRESP3  actually  calculated  absolute  count  rates,  in  terns 
of  In  i . The  absolute  rates  were  then  normalized  by  division  by  the 
predicted  rate  of  a tube  at  (0°,0°). 


* «*—  ■**#«  «-»*  *Vj>*  a*o> 
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TABLE  4 

Measured  ZAD  Counts,  Normalized  totes,  and 


Unbiased  Estimates  of  Error* 

ZAD 

N 

Reciprocity  r<<;) 

Estimated 

Position 

tow  Data,  no. 

Normalized  Rate,  KT 

Error,  i 

10®  N 

m 

1013 

.971 

* .059 

894 

1185 

1112 

734 

10°  S 

1185 

9 

714 

.982 

♦ .063 

20°  N 

1071 

571  831 

f 

.804 

♦ .051 

1036 

740  1276 

30°  N 

534 

894 

9 

301 

595 

.550 

+ .046 

30°  E 

292 

333 

65* 

9 

612 

.505 

+ .047 

30°  S 

698 

208 

1556 

9 

356 

.565 

♦ .051 

40°  N 

48 

35 

395 

9 

356 

.109 

+ .021 

40*  ME 

178 

200 

7!* 

9 

599 

.289 

+ .031 

169 

155 

40“  SE 

612 

9 

610 

.265 

+ .030 

40®  S 

62 

50 

595 

9 

624 

.091 

+ .015 

•None  of  the  corrections  discussed  in  Sections  IV  D through  IV  G have 
been  applied  to  either  the  rates  or  the  error  estimates. 
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APPFNDIX  A 

Muon  Telescope  and  ZAP  Circuitry 

The  folkwing  three  pages  contain  the  block  (i i aqrams  ■;!  the 
('iKtromc  components  of  the  wide-anqk  muon  telescope  ind  Mie 
/.on  i til-angle  dete<'tors,  exclusive.  <>i  tin  c iron  try  foi  the  auto- 
matic data  output  svstem.  The  hiqh-voltaqe  supplies  are  omitted 
from  the  plate  circuitry  and  ZAP  circuitry.  Separate  hiqh-vol taqc 
supplies  are  used  foi  the  stop  tank  tub's  (2200  v or  2100  v) , top 
and  bottom  plate  tub-..  (2200  v)  , and  zenith-angle  detector  tub's  (lb'O  v)  . 


!i  * 


3H  tube 


aEsfi 


IIIBW 


'j  nsec 
cl ipline 


oinc 


300  nsec 


delay 


DIS 


3 Msec 
width 


to  scaler  & multifile 


coincidence  unit 


com 


window 


pk 


Multichan. 


Analyzer 


Circuitry  of  stop  tank,  and  of  muon  identification.  Note: 

1.  High  voltage  supply  to  stop  tank  tubes  separate 
from  that  of  top  and  bottom  plates. 

<?.  Cliplines  cause  a pulse-shortening  by  reflection, 
riving  an  effective  resolving  time  of  5 nsec  between  5"  tuoes. 

3.  Limiters  do  not  function  as  such,  remain  in  circuit 
only  to  preserve  timing  (after  1 February  1973). 

4.  300  nsec  ielay  causes  all  M«uens”  decaying  between 
0.5  and  3. * cec  to  register  a count.  The  time-in*'frai 

of  th*-  auo«  decay  curve  over  this  interval  yield#  m auto- 
matic efficieac'*  of  7*%. 
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APPENDIX  H 


Zenith-Angle  Detector  Rato- iYodielive  Programs 

•nit*  following  pages  list  tiio  predictive  computer  porgrams 
described  m Section  III  B,  .inti  contain  same  samp  1 < • runs  of  thes- 
programs.  TVio  programs  ZADRESP,  ZADRESP?  and  ZADRESP 3 represent 
a developmental  sequence,  and  are  presented  in  chronologic  1 order 
ZAEWRITF  is  used  to  feed  the  coordinates  of  the  subdivisions  of  a 
ZAD  into  a file,  named  ZADQOQRD , from  winch  they  may  iie  read  b/ 


ZADRESP 2 or  ZADRESP  3. 


1 


H2 


ZAADRESP 

I 0 

80  ! 

PRINT  "I  NPUT  N«««"  j 

i 1 

90 

INPUT  N 

12 

110 

PRINT  * I NPUT  ZE?A,XI...“| 

13 

120 

INPUT  E„A 

l 4 

130 

LET  D = 1 

15 

131 

LET  T“2 

1 $ 

140 

LET  X<D>  = «0*SlN<E/5  7.29)*COS<A/57,29) 

17 

150 

LET  Y<D)r60*SlN<E/57.29)*SlNCA/57.29> 

18 

160 

LET  Z<D)r60*C  OS  (E/57.29) 

19 

1 70 

LET  R (D)  - 60 

20 

180 

FOR  X-*7 .5  TO  7.5 

21 

190 

FOR  Y=-7,5  TO  7.5 

22 

200 

FOR  Zr-7.5  TO  7,5 

23 

201 

LET  X(T)cX4-(X<D)-X)*(33.5-Z)/(Z  (D)-Z> 

24 

202 

LET  Y(T)tY4-(Y(D)-Y)*(33.5-Z)/(Z<D>-Z) 

25 

203 

IF  X(T)  >24  THEN  260 

20 

204 

IF  X (T)<-24  THEN  260 

27 

205 

IF  Y (T) >24  THEN  260 

28 

206 

IF  Y<T)«-24  THEN  260 

29 

210 

LET  RsSQR<(X<D)-X>**2  ♦ (Y(D>-Y)"2  + (Z(D)-Z)*2> 

30 

220 

LET  U=X<D)*(X<D)-X)  ♦ Y(D)*(Y(D)  -Y)  Z(D)*(Z(DW> 

31 

230 

LET  Q = (Z(D)-Z)/R 

32 

240 

lET  P -U/( (R  *3)*R (D) ) 

53 

250 

LET  C=C+(Q*N)*P 

34 

260 

NEXT  Z 

35 

270 

NEXT  Y 

36 

280 

NEXT  X 

37 

290 

PRINT  “RATE  =*tCx"TIMES  IC0)*SIGHA." 

38 

READY 

99  9 END 

I 

I 

I 

I 

I 

I 
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/OLD  ZADRESP/RUU 

V.)  ZADRESP  12  WAY  75  23:00 

INPUT  M * . . ? 4.7 
INPUT  ZETA.XI...7  0,0 

RATE  - 1.105  TIMES  1(0 )*SIGMA . 

9.233  SEC . 26  I/O 

2) READY 

RUN 

ZADRESP  12  MAY  75  23:01 

INPUT  N...7  4.7 

INPUT  ZETA.aI...?  40,45 

KATE  : 0.327894  TIMES  I(0)*SI',MA. 

9.275  SEC.  26  I/O 
HEADY 


* * * * * 


GENERA T ION  OF  A NORMALIZED  COUNT  AT  40  DF ORFFS  NORTH- 
EAST USING  PROGRAM  "ZADRESP",  THF  ACTUAL  NORMALIZATION 
TAPES  PLACE  BY  DIVIDING  THE  TWO  ABSOLUTE  RAT FS , NOT" 
THE  DIFFERENCE  OF  A CONSTANT  BETWEEN  TKFSF  ABSOLUTE 
RATES  AND  THOSE  CF  "ZADRESP3M.  THF  "ZADRFSPJ"  VALUES 
APE  CORRECT. 


34 


Zf'it".  >/’ 


1 ~u  i_c.T  0 * 1 
Id*.-  it.Z  I-c 

' 1 b7  Lr.T 

I 90  r'lirj 

I i SS  Phi  Ml  I .MPJT  EXPONENT  0^  COS*.  THE!  A;  . . : 

| 197  I MpjT  N 

E0O  Foh  L- 1 TO  16 
E10  hr.AU  #1:  X'D) 

-ZO  rOn  x*  ■ 7 To  7 STEP  i. 

L5u  Fv)lv  Y - *■ 8 TC  b oTlP  l. 

E40  FOh  Z--b  TO  b bTc.P  c 

| *-41  LET  X(T)-*X  t CX(D)-X>*<55.1 -Z)/(Z(D)-  Z) 

* L4i  LET  Y(T): Y+mO)-YJ*(55.W)/(Z(D)-7) 

<.45  IF  AvT)*-L5  ThLN  50C 
f c4'<  IF  X(T)>ES  T/-.EN  500 

| L46  iF  Y(T)«-L4  THEN  cOl 

E47  IF  Y(T;»t.'i  1/ic.N  50C 

L'  v l.-T  rt.‘  J<*hUXU>)-X>  *L  * CY(O)-Y)  \ ■*  i Z(  IO-Z) 

I ^0  LtT  U:X(D)»vX(''>-XMYiD>*(Y(D)-Y)^7.  * P>»  < ZCUJ-Z) 

w'v  i».'T  yi"(ZC»J)*‘Z)/i\ 

EEC  LET  P.'U/(vK**5)*h(D)) 

I <.90  LcT  C:t..4  16 

50:.  NEXT  Z 
510  NEXT  Y 

I5L0  NEXT  X 
550  NEXT  L 
540  Pit  1 NT  C 
SS9c.N0 


I 


85 


i'JTh' l ThF.  PROCEDURE  FOR  CRTA  1^1  Ml.  NPRVARI7"!>  COUNTS 
FRO,-  "ZADREOPp"  AND  "ZAORFOPr  A HP  I )F  NT  I CAL . 
Thtlitf-CKE,  OF  PA RATE  SAMPLE  PUN  OF  ”7ADR-'SPc 

AND  **ZA DWfilTE  ” 10  NOT  INCLUDED.  " ZADPEOPc" 
USES  TKE  SAMF  DETfcT T CR  tfFABl:RF|y’ENTO  AO  DOEO 
"ZAUKtSPi",  BUT  DOEO  NOT  PROVIDE  AO  VIICR 
I NFCRXAT I C.N  . 
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/OLD  ZADVRITE/RUN/OLD  7.ADRESP3/RUN 

2)  ZADWR1TE  12  MAY  75  22:5c 

INPUT  DETECTOR  ZENlTh,  AZI.-UTh  ANGLES...?  0,0 


4)  ZAD  RE  SP  3 12  MAY  7 5 22:32 

INPUT  NORMALIZING  FACTOR?  1 

INPUT  DETECTOR  ZENITH,  AZIMUTH...?  0,0 

INPUT  EXPONENT  OF  ZENITH-ANGLE  DIST...?  4.7 

AVERAGE  RAWGt  TRAVERSED  = 8.42407  INCHES. 

COUNT  RATE  = 25.0236  TIMES  I(0)*SIGMA. 

25.0236  = NORMALIZED  RATE. 

27.976  SEC.  60  I/O 
4)  READY 

/OLD  ZADWRI TE/RUN/OLD  ZAD HE SP 3/RUN 

2)  ZADWRI TE  12  MAY  75  22:37 

INPUT  DETECTOR  ZENITH,  AZIMUTH  ANGLES,..?  40,45 


4)  ZAD RE SP 3 12  MAY  75  22:37 

INPUT  NORMALIZING  FACTOR?  25.0236 

INPUT  DETECTOR  ZENITH,  AZIMUTH...?  40,45 

INPUT  EXPONENT  OF  ZENITH- ANGLE  DIST...?  h.7 

AVERAGE  RANGE  TRAVERSED  = 10.7775  INCHES. 
COUNT  RATE  = 7.03727  TIMES  1(0)*SIGMA. 
.281225  : NORMALIZED  RATE. 

cl  .504  SEC.  60  I/O 
4J  READY 


4c  * * * * 


GENERATION  CF  A NORMALIZED  COUNT  mT  43  DKG! 
HAST  COIN!  PROGRAM  ZADRESP3" , 


TS  NCYTP- 


ZADYnlTE 


1 FILE  Hi  ZADCGGnD* 

2 SCkATCH  #1 
5 LiT  DU 

10  LET  RCg): 63.0 

If'  Phi  NT  INPUT  DETECTOR  ZtNITH  AZIMUTH  ANGLES 

20  INPUT  E 

25  IF  E- C Thc.N  LOO 

30  LET  E=E/57.29 

40  LET  A=A/57.29 

50  LET  XO  "h  (D  )*  SI  N(  E) 

51  LET  XQ-XO+C*.‘> 

80  LET  Dl=l.l*COS(E> 

85  LET  02  -1.1 


loG 

FOh  Xl  = X0-t.5*Dl  TO  XO+l. 

5*D1  STEF 

140 

FOR  Yl=- 1.5*02  TO  l.5*D2 

STEP  D2 

150 

LET  Z:S(*R(k(D)*2  - XI  *2  - 

Y1  *2) 

152 

LET  W=aTN(Y1/X1> 

154 

LET  Rl  = S<ih<Xl  *2+Yl  *2) 

156 

LET  Y:hl*SlN(W+A> 

1 5fc 

LET  X-R1*C0S(W+A) 

160 

toftITL  flsx 

170 

NEXT  Y l 

160 

NEXT  XI 

190 

STOP 

200 

LETZ-  63.0 

210 

FOh  X--1.15  TO  2.15  STEP 

1.1 

L2Q 

FOR  Y-* 1 .65  TO  i .65  STEP 

1.1 

230 

IkRITE  MsX 

240 

NEXT  Y 

250 

NEXT  X 

999 

END 

89 


APPENDIX  C 

Muon  Telescope  Analysis  Programs 

Program  MLJRESP  is  a direct  offshoot  of  the  ZADRESP3  program. 

Hie  only  change  in  detector  description  is  the  replacement  of  the 
ZAD  area  and  solid  angle  with  the  subdivisions  of  the  top  plate.  A 
factor  of  cos  is  required  due  to  the  top  plate's  oblique  angle  with 
respect  to  the  incoming  flux  away  from  the  vertical 

Programs  DETECTOR,  whose  detector  dimensions  are  only  approximate, 
and  DCTBCTR2,  where  the  dimensions  are  considered  accurate  to  ±0.5", 
are  the  products  of  the  "trajectory"  method  of  predicting  detector 
response.  They  are  not  suitable  for  verification  by  zenith-angle 
detectors,  as  is  MURESP.  However,  they  employ  the  same  detector 
description  as  does  MURESP;  further,  multiplying  the  curve  in 
Figure  15  by  I (8,4.) sin  8 and  integrating  over  the  wide-angle  telescope 
aperture  should  give  the  same  results  that  MURESP  gives. 


MUrti-SP 


10 

1' ft  1 U T INI  JT  u> 

'POfltMT 

Of  l8Ct#  I‘.- 

K> 

IMPuT  »» 

17 

PrtlNT  JCJT.l  ;.ALKU) 

Oh 

No?^7i1  HM*.  » w 

*.0 

INPUT  ,t 

IF  h-L  The*  cl 

« T 

LcT  As-LQ.V 

iu4 

LET  fc-'cS: 

Li> 

•iOTO  aO 

27 

LET  A:3.5 

LB 

LET  a- d #5 

30 

LeT  l! 

<»G 

FOn  XI  =«  T J ti 

TEP  t 

» 

!>0 

FOR  Y1 -“Ll  TO  i 

c\  STcP 

6 

60 

FOR  XL  = -7  TO  7 

STEP 

* 

c 

70 

FOR  YL2-8  TO  8 

STEP 

(L 

SO 

FOR  Z2  = -8  TO  8 

STEP 

C 

50 

LET  RrSLhUXL- 

XI ) *2 

+ 

iYL-YI  > *'c  ♦ 

100  LET  C=(Z1-Z2)/R 

110  LET  jr&*G6*(C“lN*l  ))/{h\) 


!LC 

LET  T:T+0 

130 

IF  C<COS(l 0/57.29) 

TmEN 

16C 

l AO 

lET  F(1):F(1)+D 

1 50 

UOTO  300 

19C 

160 

IF  C<C0SIL0/57.L9) 

T heft 

17C 

LET  F(L)2FCL)+D 

180 

GOTO  300 

L.iu  0 

190 

IF  C«COS(30/57.L9> 

T hi.  M 

200 

LET  F(3) : F(3}+  D 

c 1 0 

goto  oOO 

L2u 

IF  C«COSlAO/57.L9) 

THefc 

i.3v 

LLT  F!A)  = F(A)0 

LAO 

uOTO  3G0 

L5o 

LET  FO):F-.‘3>+D 

ECO 

LET  t»l  ••  oi  1 ♦ 1 / 

31U 

let  S=G*1 

3L0 

next  zl 

330 

NtXT  YL 

C>*G 

NEXT  he 

3?v 

NEXT  Yi 

3bC 

NeXT  XI 

PkUT 

iby-y  < 
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tURESP  (CONTINUED) 


370  PrtlNT  “ TOTAL  COU'.T  -.aTE  - “ ; I;"  T Ir.Eo  I(C)*oMP,A.' 
3tJu  PitlUT 

5'o‘j  PRINT  LOWER  oOUND  f JPPcK  uCuOD  t CuJnTO 

330  FOK  l:l  TO  5 

400  PRINT  10*a-i)>10»L,F(L) 

4i  G NEAT  L 
4L0  PhlNT 

4iG  PRINT  "AVERAGE  jLCANT 
999 iND 
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i~C  Fuh  Al»  • “O  1 j 6*j  oT— P j 
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81  Li. 7 Y^Y^  * L/ c •♦  (9't- 
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l Jv-  i ► \><.A  T v.  v j 

110  IF  Y«-Za  T.W.N  4.00 
U0  IF  Y>c.ii  TriEN  ZOO 

r. o let 

cCiO  Nl/J  Z' 
i.  SfJ  NEAT  Y? 

Z4G  NEAT  Ar- 
2.50  NEaT  P 

7.00  PnlNT  A7*i>*ICG/c;i6/.»)  *3) 
539  ENL) 


oTEP  L 
ET-P  ., 
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4W>3/2>*  TAN(T/‘  7,,.0>Cu. 
(9*«*  ZL-S/Z:»TAWa/L7. J I! 
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DE  TtOl«£. 


U ‘-i-.lCH  ZLMlTh  ANOLe" 

..C  Itil  Jl  £ 

Hn  A--C  TO  ^40  OTEP  EC 

Au  rOh  XE---7  TO  7 STEP  E 

>c  F Oh  Yi.r-8  TO  b STEP  e 

6u  rjh  Zk:-6  TC  8 STEP  L 

TO  Lil  XI  - XZ.+  CiE.tS-ZE) 

'60  LET  Y 1 = YE*  ( 5e  *Ef>  • Z?  ) 

SO  If  Xl<-k;5.i.  ThEN  140 
100  IF  XI  >04.5  T.-ttN  1 ; j 
i 1 v lr  Yl<-E4  f /»£.  1W  140 
*cii  IF  7l>Lu  TnEN  140 
1<50  LET  C' v-n 
i'ic  .Mr.  XI  U. 
it>0  NeXT  Y i. 

IU  NEXT  XL 
1 7 nEXT  A 
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APPBDIX  D 

Statistical  Analysis  Program 

Program  RESULTS,  listed  in  the  following  pages,  has  two  output 

forms.  First,  it  may  be  used  to  print  out  a table  of  reciprocity- 

normalized  data  points,  corrected  or  uncorrectnd , along  with  the 

estimated  errors  and  the  raw  rmnbers  of  counts.  Second,  it  can 

evaluati  the  flat  and  weighted  S values  for  the  predicted  data  listed 

in  lines  400-410,  for  either  corrected  or  uncorrected  data.  The  S 

values  aenerated  by  excluding  the  points  at  40°  N and  40°  S may  be 

obtained  oy  the  following  modification: 

320  DATA  8 

560 

590 

410  DATA  (predictions  for  40°  NE,  40°  SE  only) 


9'i 


RESULTS 


1 *•  I t't  , P k-0 ) j Pf  # ■ *.  w / 

*.  ;-mI  fc.z.‘. 

3 MAT  F:  *.U. 

A PhlNT  COhnUlLU  1)  0>i  OtiCG.-thECIEUl  ' ; 

> I MPdT  <* 

0 HEhD  \ 

*3  hLAl/  L,Y,Z 

1 6 IF  L-i  Thtrt  U 

17  IF  (,:L  THEN  It) 

18  LET  E< L ) - Y 

19  UT  F(L>-Z 
U GOTO  It/ 

il3  FOP  L - i T : V 
30  he  AO  s'  \ i.  .* 

:•>  *o:  u 

AG  FOh  L-l  TO  N 

43  IF  L-o  THEM  1100 

SO  HEAD  A<L>,ba>,CU..\D<L> 

M LIT  A.-m’.U 
3i.  U T L ’ *>  * 1. ) 

3 i U T C - 1. 1 L ■ 

34  LET  D U L ; 

80  LLT 

83  LeT  hi  L)  "n  »L)*  Z(L  > 

90  LET  hi  - Si* H 1 l * i > U.-. . A ' • I /o*h  > u /'•♦  1 / .'Uhl  C)  + 1 /L>L!t v j / ~ ’■ 
IOC  i.  *,  T il! ' I * Ouh  « 1 ■*  i / -j  v»i<  l A ) “ 1 . Li*  hi  L ) “ 1 / OQh  ( C / “ « # jx’  n<  < < 

1 10  LET  &iL'  fi(  l.)»  El  ‘rU) 

uo  LiT  y-Wf  (h(L)-P(L>) 

130  LET  Cl  UI  + SUUdIDOCL^M  hll  >-i  ID)  UT.L  ) 

1 AO  MEAT  L 

139  PHI  NT  “DO  YC.J  W h;»:  I a oL  E i 1 ' oh  ' ; 

100  INPUT  0 

i 7J  IF  UL  T/iFM  3 CO 

18C  PniNT,,  hcflPhOCITY 

JdL  PriUT  "h#*W  DATA/,"  M/nf-AL  IZ  .0  _ 

18A  Phi  NT  J’OailoCIUN.u/^AT:-  , LlGr-'A  „ 

U6  Phi  NT  * ",  " ■*•  

ISC  PGh  L-l  TO  ■« 

LOO  hc.AD  LI  ^ 

Eli  PnINT.A vL);  / }ciU/ 

*.13  if  l-3  Then  ;l3o 
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